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Abstract 

For neutrinos of ©(lOkeV) energies, their oscillation lengths are less than a few 
hundred meters, thereby suggesting the fascinating idea of oscillation experiments 
of small geometrical size. To help evaluating this idea, a formalism is developed for 
calculating the neutrino ionization cross sections for H as well as the noble atoms. 
This formalism is based on the use of spin-independent atomic wave functions, and 
should very accurately describe the ionization spectra for H, He, Ne and Ar. The 
accuracy is considerably reduced for the Xe case though, where the spin dependence 
in the wave functions is non- negligible. Nevertheless, even for Xe the results remain 
qualitatively correct. In all cases, the atomic ionizations cross section per electron 
is found to be smaller than the neutrino cross section off free electrons, approaching 
it from below as the energy increases to the 100 keV region. At the 10-20 keV range 
though, the atomic binding effects in the cross sections and the spectra are very 
important and increasing with the atomic number. They are cancelling out though, 
when total ionization cross section ratios, like i^^/i^e or v^/ve^ are considered. 



PACS numbers:13.15.+g 



1 Introduction 



Neutrino interactions have been observed and analyzed at energies from a few MeV (in 
reactor and solar neutrino experiments) [Tj, to several hundred GeV at accelerators. The 
lower keV energy range though, has not attracted much interest up to now, due to the 
smaller cross sections there. 

Nevertheless, when intense neutrino fluxes from radioactive decays become available, 
the possibility of very interesting experiments may arise, such as the testing of the solar 
neutrino oscillations using terrestrial experiments of rather small geometrical size [2]. To 
see this, we remark that recent results from the KamLAND experiment gives 

Am?2 = 8.2toi X 10~^ eV^ , tan^ ^12 = O-Sdtofr , (1) 

which would suggest a corresponding neutrino oscillation length of about 15 meters for 1 
keV neutrinos, increasing to 150 meters when the energy reaches the 10 keV level. The 
oscillation length associated to Am^g ~ Am^g is about 27 times smaller 

To perform oscillation measurements at keV energies, the ionization of atoms by neu- 
trinos may be used, taking of course into account the binding of the electrons. As an 
example we mention that the realization of an intense Tritium source may allow the 
study of z/g oscillations in a room size experiment. 

Such experiments may also help in investigating the weak interaction couplings at very 
low energies |2j, the measuring of sin^ ^13 (5j, the search for measurable contributions from 
a neutrino magnetic moment |21 El El , or contributions from possible new gauge bosons 
jS]. These are interesting questions, which were never investigated at very low energies, 
and may provide useful additional information on neutrino physics. 

On the theoretical side, there are three calculations known to us investigating such 
effects. In [H], the energy spectrum of electrons knocked out from ^^F{Z = 9) and 
^^Mo{Z = 42) atoms is computed, using relativistic electron wave functions. It is found 
that the electron ionization spectra differ significantly and are always smaller than those 
of the free scattering case. 

In jH], we have also computed the total cross section on free electrons, as well as on 
bound electrons in some light atoms. At that time we studied the three atoms H, He 
and Ne. The result again was that the ionization cross section per electron is always 
smaller than the cross section from free electrons [Hj, and never larger Moreover, this 
ionization cross section always decreases as we proceed from H to He to Ne, obviously 
due to the increase of the binding [U]. 

For planning and carrying out these new experiments though, it is needed to have the 
energy-distribution of the ionization electrons for the noble gases He, Ne and possibly Ar. 
The reason is that the detector sensitivity may generally depend on the ionization electron 
energy and it anyway diminishes below a certain limit; say e.g. ~ lOeV 0. In principle, 
Xe could also be useful, but such an experiment would be much more expensive than a 
Ne or Ar one [3]. Finally, the use of Kr is hindered, by the existence of a /^"-emitting 
isotope 
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The aim of the present paper is to present a method for calculating such spectra. 
We restrict to cases where it is adequate to consider the very accurate spin independent 
atomic wave functions of ^I]. This way, the most interesting cases of the He, Ne and Ar 
atoms are covered, for which the results should be quite precise. 

In principle our method could be extended to Xe, for which spin independent wave 
functions are also given in ^T]. But the results would then be less accurate^. This 
is suggested by the observation of a considerable spin dependence in the measured Xe 
binding energies indicating that a dedicated relativistic treatment is needed for it. 
Thus only a very brief discussion of the integrated Xe ionization cross section is offered 
here. 

Finally, since the properties of the detectors needed to experimentally study these 
ionization spectra are intensively investigated at present |21EIII21, the feasibility of this 
task is conveniently left outside the scope of the present work. 

At the very low energies we are interested in, the Standard Model dynamics described 
by the diagrams in FigQ induce the effective interaction Lagrangian 



Gp 

'71 



_ , (1-75) 

T^eT ^ T^e 



+ 



- 0(1-75) , - p(l -75 



(2) 



describing the electron interactions with any neutrino flavor. The Ue and {up, Vr) couplings 
are respectively given by 

Ve = l+ , ae = 1 , 

Ve = -1+ , fie = -1 , (3) 

while Gi? is the usual Fermi coupling. 

We first concentrate on the Ve{i^e) cases. From Q, the squared invariant amplitude 
summed over all initial and final electron spin-states for the process 

Ue{Pi)e-{P2)^Ue{P;)e-{P^) (4) 

is calculated, with the four-momenta being indicated in parentheses. The various particle 
energies are denoted below by Ej, and the standard variables 

s = {Pi + P2f>ml , t = {P^-P;f , u = {Pi-P^f , (5) 

are used. We note that (H)) describes the Pe scattering from either a free or bound elec- 
tron, the difference being determined by £"2; and, of course, the folded in momentum- wave 
function in the bound electron case. 



^For an estimate of their error, see the discussion at the end of Section 2. 
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Neglecting neutrino masses, we start from the case where the initial electron is free, 



so that Pf 



p2 
^3 



0, P| 



p2 
-'4 



with me being the electron mass. Summing over 



all initial and final electron spin states, we then have 
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(6) 



Due to crossing, the corresponding \F\^ expression for the Ue case is simply obtained from 
the Ve result, by interchanging s m in 

In the lab system where the initial electron is at rest {E2 = me), the differential cross 
section describing the energy distribution of the final electron is^ 



dcr(z/ee 



free 



'^^'^F J + aefE^ + {Ve - aef{Ei + me - E^f 



SnEf 

+ me{vl - al){me - E^)^ 



(7) 



which, after it is integrated over the allowed range 



rrie < E4 < rrie + 



2Ef 



nie + 2Ei 



leads to the Ue total cross section off free electrons 
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(9) 



Starting from ((7j), we present in^ FigOthe energy distribution of the final electrons, 
in Ue and scattering off free initial e~. In these figures Ee = E4 is the final electron 
energy, and E^, (Ep) denote the incoming neutrino (antineutrino) energy Ei. As seen 
there, the energy distribution always has maximum at small Ee, which is also seen in all 
z/e(z/e)-induced ionization cases; see below. 

In Section 2 the formalism describing the ionization of atoms through neutrino (an- 
tineutrino) scattering is presented. This is done first for the Hydrogen atom, and it 
is subsequently generalized to any atom characterized by complete electronic shells. A 
discussion of the accuracy of our method for the He, Ne, Ar and Xe atoms is also given. 

In Section 3, we present the results for the energy distributions of the knocked out 
ionization electron in Uei^^e) scattering for H, as well as for He, and Ne, using the analytic 
spin- independent Roothann-Hartree-Fock wave functions published in |TT]. In addition. 



^Notice that the laboratory angle of the emitted electron is not an independent variable, but is fixed 
completely by its energy. 
^We always use = 0.23. 
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we also present results for Ne ionization through or v^. scattering, which might be 
generated through z/g oscillation. Results for the integrated cross sections are also given 
in the same Section 3, where we also include our predictions for the Xe case. The Con- 
clusions are presented in Section 4. Important details on the kinematics and atomic wave 
functions are relegated to the Appendices A.l and A. 2 respectively. 



2 The formalism. 

The ionization of atoms through neutrino scattering off atomic electrons, requires from 
atomic physics the electron binding energy and the wave functions in momentum space. 
The energy of the initial bound electron is fixed by the binding energy e, so that 

E2 = me + e , (10) 

with e being negative. For the simplest H atom, e is fixed by the Balmer formula, while 
for more complicated atoms, we use the experimental measurements ^j. 

The three-dimensional momentum of the bound electron however, (compare definitions 
in ()A.1|) ). varies according to the probability distribution 

hT, ( a ^ M2 Pidp2dcosg2d02 . . 

|^nim(P2,^2,02)| , (H) 

determined by the momentum wave function'^ ^E'n«m(P2, ^2, 02), when electron spin effects 
are neglected. Therefore, the bound electron may get off-shell with an effective squared- 
mass given by 

m^ = P^ = El-pl , (12) 
so that, the standard kinematical variables defined in (0) satisfy 

s + t + u = fn^ + ml . (13) 

The squared invariant amplitude for the i/g-electron subprocess, summed over all initial 
and final electron spin states^, is written as [H] 

+ {Ve - aef{u - ml){u - m^) + 2ml{vl - a^)t| , (14) 

which of course becomes identical to (jHl) when w? — > rp?^. As in the free electron case, the 
antineutrino result may be obtained from ()14|1 by interchanging s *^ u. 



^The momentum here is connected to the wave function determining its probabihty distribution. It is 
the conjugate variable to the position coordinates, and their wave functions are related through a Fourier 
transformation. See Appendix A. 2. 

^Justified because we neglect any spin dependence on the wave function. 
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For clarity, we first consider the ionization of a Hydrogen atom being initially in its 
ground state. The energy spectrum of the ionized electron is obtained by averaging over 
the bound electron momenta according to its wave function (see Eqs.(15) of 0), and 
subsequently changing the dcr/dw-distribution to da/dE^, which brings in the derivative 
du/dE^. We thus get 

dc^H _ 1 f d02d cos 02 Pldp2 I . ^|2|p|2dM 

d^ - 64^^ j (2vr)3(. - m^) l'^^^"^^^^' dE, ' ^^^^ 

The variables {p2, O2, ^2) are the momentum and angles of the bound electron in the rest 
frame of the atom, (compare ()A.1|) ): and^ ^100(^2) is the ground state momentum wave 
function defined in ()A.16|) . Finally E2 is determined through ()10|) , which for the H ground 
state (n = 1, / = 0) is given by 



e 



= _::!£l^ = _13.6eV^ . (16) 



The kinematics are fully explained in Appendix A.l; (compare (^). According to it, 
for any (^2,02) in the range ()A.3|) . and any value of the bound electron's momentum p2 
and the incoming neutrino energy Ei, the polar angle of the ionization electron ^4 is a 
function of its energy i?4 given by^ flA.6|) . This function is used to determine through 
()A.9|) . the expression for dM/d£'4 needed in (fT^ 

We also note that the angular dependence of the integrant in (fT3jl is only due to ()A.9|) 
and the effect of ()A.4|) on The angular integration is done numerically, with its 

range fixed by ()A.3|) . For the numerical evaluation of the p2-integral, the relevant part 
of the p2-range is determined by the form of the electron wave function, as discussed 
immediately after ()A.16|) and at the end of Appendix A. 2. 

We next turn to the general case of any of the noble gas atoms He, Ne, Ar, Kr and 
Xe. Their wave functions are discussed in Appendix A. 2 JTj. Since for noble atoms all 
electronic shells are complete, a summation of the form 



_ 2/ + 1 



|2 



always appears in the wave function contribution to ionization, washing out any angular 
dependence from it. As a result, the energy distribution of the ionization electron for any 
noble atom, normalized to one electron per unit volume, may be obtained from (fTK|l by 
replacing 

|^ioo(p2)r ^ -^^!.2[R,0{P2)? + 2[R20{P2)? + 6[R21{P2)? 



^This wave function has no angular dependence, because of the vanishing of the orbital angular 
momentum. 



^See the discussion immediately after ljA.6|l for resolving any ambiguity. 
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+ 2[^30(P2)]' + 6[i?3l(P2)]' 

+ 10[^32b2)]' + 2[i?40(P2)]' + 6[i?4l(p2)]' 

+ 10[RM]' + 2[i?50(P2)]' + 6[i?5l(P2)]'| , (17) 

where Rni{p2) are the radial momentum wave functions defined in ()A.13|) . and Z is the 
atomic number. 

Focusing to the expression within the curly brackets in the r.h.s of (fTTj) . we remark that 
by restricting to just the first (one, three, five, eight, all) terms, we obtain respectively 
the results for the (He, Ne, Ar, Kr, Xe) atoms, provided the appropriate Z value is used. 
The corresponding atomic wave functions are discussed in Appendix A. 2 and [TT] . 

We next turn to the discussion of the reliability of our calculation. We first recall that 
a very high accuracy of more than 8 digits has been claimed by the authors of [TTj, for 
all RHF (spin averaged) energies of their results, for all atoms from^ He to Xe. The main 
weakness comes therefore from the neglect of spin effects in the atomic wave functions. 
This should be no problem for the lighter atoms up to Kr, for which the experimental 
data on the electron binding energies indicate that the spin effects are smaller than ~ 1% 
[in]. Therefore, our results for He and Ne, and those that could be obtained by applying 
our procedure to Ar and Kr, should be accurate at the 2%-level, their main error coming 
from electromagnetic radiative corrections and small relativistic spin effects. 

For Xe though, the spin dependence of the electron binding energies shown in Table 
1 is more significant, suggesting that the error in the individual ionization cross section 
could possibly reach the 10% level^. In the discussion of Section 3 we find indications 
though, that this error is in fact diminished in the ratio of the v^/ total cross sections 
determining 9i^. The electron binding effect seems cancelling in this ratio. 

Table 1: Relative spin effect in the binding energies 
of the various electronic quantum states in Xe, P^ . 



Quantum states 


Spin effect 


2pi/2 - 2p3/2 


6.6% 


3pi/2 - 3p3/2 


6.4% 


3(^3/2 - 3(^5/2 


1.8% 


4pi/2 - 4p3/2 


0.8% 


4^3/2 - 4^5/2 


2.9% 


5pi/2 - 5p3/2 


10.4% 



^See in particular the last sentences in Appendix A. 2 quoted from their paper. 

^Note that the stronger singularity of the relativistic wave functions at small distances, which increases 
the momentum wave functions at large momenta, should have no effect on the ionizations, since in this 
range the wave functions are strongly suppressed. 
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3 Results for neutrino ionization of Atoms. 

Using (|T31 ll6 |) and the H-wave function in ()A.16|) . we obtain the results in Fig 121 describing 
the energy distributions of the knocked out electron in Uei^^e) ionization of Hydrogen, nor- 
malized to one electron per unit volume. The results apply to various incoming neutrino 
energies. 

The corresponding distributions for He and Ne atoms are shown in FigsE] and re- 
spectively, using the wave functions of Appendix A. 2 and the parameters tabulated in 
For the electron binding energies in the various atoms (including those of Xe we 
discuss below), we use the experimental values [12] 



He : 


.He _ 
^Is — 


-24.6e1/ , 










(18) 


Ne : 


Ne 
^Is — 


-870.2e1/ 


Ne 

1 ^2s — 


-48.5e1/ 


Ne 
1 ^2p ~ 


-21.7eV 


(19) 


Xe : 


.Xe _ 
^Is — 


-34561e1/ 


.Xe _ 


-5453e\/ 


Xe 


= -A893eV , 






.Xe _ 
<^3s — 


-llAQeV , 


.Xe _ 


-961eK , 


.Xe _ 
^3d — 


-681AeV , 






.Xe _ 


-213.2eV" 


.Xe _ 


-Ub.leV 


Xe 


= -68.3eV" , 






.Xe _ 


-23.3ey , 


.Xe _ _ 
^Sp ~ 


-12.5e1/ . 






(20) 



As seen from FigsI2]Elthe energy distribution in all z/e(z/e) -induced cross sections is maxi- 
mal at small Ee = E^. 

Concerning FigslHHHElwe may remark that for Ei in the few lOkeV range, the differen- 
tial dcr/di?4 ionization cross sections per electron are always at the level of l{)~^^cm?eV~^ . 
Moreover, they are always smaller than the "free electron" ones, approaching them from 
bellow as Ei increases; see the (c) and (d) parts of each of these figures. For a fixed 
E'l-value though, dcr/di?4 always decreases as we go from Hydrogen to Ne and then to 
Xenon. 

The integrated z/e(z/e)-ionization cross sections for final electron energies ii^4 > mg + 
lOeV^, are shown in FiglHl Here again the cross sections are normalized to one electron 
per unit volume. In the same figure the corresponding integrated cross section off free 
electrons are also included, as well as our predictions for Xe. As seen from FigE^i-jb, all 
cross sections are of the order of 10"^"^ cm^ at Ei ~ 15keV, increasing with the neutrino 
energy. The ratios of the integrated ionization cross sections to the free electron one are 
shown in Figs|nt,d. As expected, the integrated ionization cross section is also always 
smaller than the free one, approaching it from below. For H and He the approach is very 
fast, but it becomes much slower for heavier atoms; see e.g. the results for Ne and Xe. 
Thus at El ~ 50keV the Ne cross section is still about 5% smaller than the free e-one, 
while for Xe the decrease is at the ~ 30% level. 

Such integrated ionization cross sections for H, He and Ne have already appeared in 
p. The present results are consistent with those. We should point out though that in 
jU], the integrated cross sections were calculated over the entire physical region E4 > mg. 
Moreover, in that work we had used a very rough approximation for the He wave function 
based on the Zeg-idea, while the Ne results were based on the old fit of Tubis ^1]. Thus 
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although the present results are consistent with those of 0, prospective users should rely 
more on the present ones based on TT. After all, it is only here that the energy distri- 
bution of the ionization electron appears. 

We have already estimated that for Ve oi in the 10 keV level, the oscillation length 
is about 150 meters; see (0). Consequently, the other neutrino flavors should also be 
generated, as the initial z/g or z/g proceeds in a volume filed with a noble gas. Thus, in a 
Tritium decay experiment, and should appear in relative amounts determined by 
the mixing angles and the distance from the source. 

The treatment of the effects is exactly the same as for z/g, the only difference being 
that in (0 inj) we now have to use the second set of the neutrino vector and axial couplings 
listed in (jHl); compare 0. The resulting differential cross section for scattering off a free 
electron is then given FigEJ which is strikingly different from the corresponding vJ(Vf^ 
result shown in Fig|21 A similar situation arises also for Ne ionization; compare the Vy^ 
induced ionization shown in FiglHl to the corresponding v^. effect in FigjSl Again, the E^. 
distributions tend to have a local minimum at low final electron energies in Figs l7l8l in 
contrast to the local maximum if Figs|21El 

We have moreover found that for neutrinos (antineutrinos) at the 10-20 keV range, the 
respective ratio of their integrated cross sections^° z^^/t'e (^/x/^'e), remains almost constant 
and equal to 0.42 (0.44), as we move from the free electron case, to the Ne and the 
Xe ionization ones. This suggests that all binding effects in the individual v^. and z/^ 
integrated cross sections, are cancelling in their ratio; for which the simple expression Q 
is adequate. 

It is worth remarking also on the basis of Figs|7| and |H1 that the z/^ and z/^ results 
are almost identical. This is due to the fact that ?)e — 0, for the value of the Weinberg 
angle we use (compare ©), which makes the squared amplitudes in (01111) almost s m 
symmetric. It seems therefore that a Tritium experiment may also give information on 
the value of the Weinberg angle at the keV scale ^2] ■ The results are of course identical 
to those for compare ©. 

We next briefly address the problem of the flavor oscillation of the neutrino induced 
noble gas ionization. Since for reasonable gas-densities the vacuum oscillation treatment 
should be adequate, the oscillation probabilities may be written as 

P(z/g -> z/^ + z/.) = sin2(2^i2) cos^(^i3) sin^ + sin2(2^i3) sin^ (^) , (21) 

P(z/e^Z/e) = \-P{y^-^v^^v;) , (22) 

assuming three active and no sterile neutrinos. Here 6'i2, 6*23? ^'13 are the 3 mixing angles, 
while Ajj = m| — denote the mass differences between the neutrino masses satisfying 
lAial ~ IA23I 3> IA12I [Ullj, and L is the distance from the source. We note that there 
is no dependence on the neutrino CP-violating phase 5 in (^1] 1^ . and that these same 
formulae describe antineutrino oscillations also. 

-^"Integrated for — rrie > lOeV. 
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Thus, as an initially produced z/g transverses e.g. a Ne target, the ionization cross 
section per electron varies with L as 



dcr 



Ne 



dE. 



[1 - P{ye 



V, 



(23) 



in the r.h.s of which (|151 fTTj) should be used for v = v^, and v = respectively. Corre- 
spondingly for antineutrinos. 

Using the experimental Aj^ values, we find that the oscillation length of the first term of 
PiVf. Vf^ + Vr) (see fl21|) ) is about 150m for an incoming neutrino energy of Ei ~ IQkeV , 
while the oscillation length of the second term is only about 5.6m. Since Oy^ is known 
to be very small, the picture created by ()23|1 will then consist of a few hundred meter 
oscillation, modulated by a much weaker one of a few meter size jlj. The strength of the 
modulation is only determined by ^13, and the ratio of the z/^ to z/g cross sections. As we 
have already stated, an enhancement of the theoretical accuracy of this ratio immediately 
implies also an increase to the ^13 sensitivity. 

Depending, therefore, on the achievable experimental accuracy, the study z/g oscil- 
lations in Tritium decay may help further constraining the neutrino mixing angles and 
masses. This would be most interesting for ^13, for which information might be derived if 
a future experiment manages to be sensitive to both oscillation lengths governing ()21|). 



4 Conclusions 

In this paper we have presented a formalism for the ionization of atoms by bombarding 
them with neutrinos of any flavor in the keV energy range. The interest in this energy 
range originates from the fact that the oscillations lengths for neutrinos in the few keV 
range become rather small, allowing the possibility of studying the oscillations observed in 
the solar neutrino and KAMLAND experiments, by means of a terrestrial experiment of 
small size. A highlight of such an experiment is to improve the constraint on ^13. Neutrinos 
in this energy may be obtained from various possible beta decays; most notably Tritium 
decay producing antineutrinos. 

Motivated by this, we have undertaken the present extensive study of the ionization of 
H and the noble atoms by neutrinos. To this purpose, we have developed a method based 
on spin independent atomic wave functions, for which we use the very accurate Roothann- 
Hartee-Fock (RHF) wave functions listed in ^T]. The energy of the bound electron in 
each atomic state is fixed by its binding energy, while its momentum varies according to 
the distribution determined by the momentum wave function, thereby generally forcing 
the electron to get off-shell. The method is very accurate for treating H and the most 
interesting noble gases He, Ne and Ar. 

On the other hand, it is less accurate for Xe, where considerable spin dependence 
appears. Thus, for Xe we have simply given the results of the present method for the 
integrated ionization cross section, whose accuracy should be lying at the 10% level or so. 



9 



This completes the overall view of the neutrino ionization of the noble gases. If the use 
of Xe in such an experiment is finally decided though, than of course a special treatment 
would be needed, based on spin depended atomic wave functions. 

We next proceed to summarize our extensive results. We have found that at the 10-20 
keV neutrino energy range, the differential cross sections dcr/dE'e are at the 10~^^ cm^/eV 
level for z/e(z/e) -induced processes, while the integrated cross section is of the order of 
10^^'^ cm^. At the 10-20 kcV neutrino energy range, atomic effects arc very important 
and cannot be ignored. They reduce the Ne ionization cross section per electron by 
almost 20%, while for Xe the reduction reaches the factor of two level. Of course, as 
the incoming neutrino energy increases beyond e.g. the 100 keV region, all these "per- 
electron" ionization cross sections approach the neutrino cross section off free electrons, 
always from below. 

We have also compared the induced reactions, with those induced by or i/^; the 
later two being equal to each other. The difference between the two comes from the fact 
that i/e reactions involve both charged and neutral current diagrams, while for f^, i/-,- 
only neutral currents contribute. As a result, the energy distribution of the final electron 
in the z/^ reactions tends to have a local minimum at low electron energies, in contrast 
to the local maximum expected for the corresponding i>e effect. Moreover for 10-20 keV 
neutrinos or antineutrinos, the ratio of their integrated cross sections is found to be largely 
independent of any atomic binding effects. 

As an overall conclusion we may state that if such small cross sections become mea- 
surable one day, neutrino atomic ionization experiments may be useful for testing the 
electroweak theory at keV energies and studying the neutrino mixing. 
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Appendix A.l Kinematics 



In the rest frame of the atom, defining the i-axis along the direction of the incoming 
neutrino, and the x2;-plane as the plane where the (final) ionization electron lies, we 
write, (compare (jH)) 



1 








-'2 



/ E2 

P2 sin 02 cos ( 
P2 sin 62 sin < 
\ P2 cos 62 



\ 



-'4 



P4 sin 64 


y P4 COS 64 j 



(A.l) 



where E\ is the energy of the initial neutrino; (i?2,P2, ^2, ^2) are the energy, momentum 
and angles of the electron bound inside the atom; and [E^^p^^Qa) are the energy, mo- 
mentum and polar angle of the freely moving final ionization electron. By definition, the 
range of these angles is 



< ^4 < TT , 

< ^2 < TT , < 02 < 27r . 
The standard variables defined in (0) become 

s = + 2^1(^2 -P2 cos ^2) , 
t = ml — 2Ei{E4 — p4 cos 64) , 

u = 2Ei{E4 — P4COS64 — E2 + P2COS62) 

where the definition (fT^ is used, and (fT^ is of course satisfied. 
The requirement of 



(A.2) 
(A.3) 



p2 



(Pi + P2 - A 







(A.4) 



(A.5) 



implied by the negligibly small neutrino mass, leads to an equation determining cos 64 in 
terms of E4 and Ei. In the free electron case, this is linear in cos6'4, and it can be solved 
immediately leading to ((Tj). 

For bound electrons, however, ()A.5|) leads to a quadratic equation in cos 6*4, which 
also depends on the momentum and spherical angles of the bound electron. It is then 
important to discriminate between the two mathematically possible solutions, among 
which only one is physically acceptable. To do this we first note that the general solution 
of ()A.5|) . may be written as 



tan,^)^^i±vS 
2) 2C3 



(A.6) 



where 



Ci = 2p4p2 sin 6'2 COS 02 , 

C2 = 4|[E4(Ei+E2)-e]'-P4(^l+P2COS^2)'} 

Cs = E4{Ei + E2)-i + P4{Ei+p2Cose2) , 



(A.7) 
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e= -^1^2 008^2 + ^1^2 + ^^^ ■ (A.8) 

A detail study indicates that whenever both 6*4 solutions of ()A.6|) satisfy ()A.2|) . the phys- 
ically acceptable one is given by the upper (lower) sign of ()A.6|) . depending on whether^^ 
Ci > (Ci < 0) respectively. This physically acceptable solution, is by its very defi- 
nition, a continuous function of the ionization electron energy, and the bound electron's 
momentum and angles entering the integration in (fT3j) . 

Once the physically acceptable ^4 solution in ()A.6|1 has been identified, dM/d£'4 enter- 
ing (jl5p is determined from 



^ d6'4 . ^ 

—- = -2^1 1 cos6'4 +P4-— sm6'4 

dii/4 V p4 d-C/4 

d^4 i?4[£'i cos6'4 + p2(cos^^2 cos 6*4 + sin6'2 sin 6^4 cos 02)] — (-f'l + -£'2)^4 ('A 9) 
d£'4 p| \E\ sin 6*4 + •p2 (cos Q2 sin 6*4 — sin 6*2 cos 6*4 cos 02)] 



Finally, the requirement Ci > C2 (note the square root in ()A.6|) ) imply that the range 
of i?4 is determined by 

Ef <E4<Ef , (A.IO) 



with 



-C/4 



i[Ei + E2) ± ii^e~ml\{Ei^E2f-(i 



(El + E2Y - i\ 



(El +p2Cos6'2)2 + P2sin^^2cos2 02 • (A-H) 



Appendix A. 2 The Atomic Wave Functions 

The atomic wave functions in the momentum space needed in ()15|) . are related to the 
coordinate wave functions by 




(A.12) 



where nlm denote the usual quantum numbers characterizing atomic states. The radial 
momentum wave function defined by the first function in the r.h.s of^^ 

^nU^) = {-l)'Rnl{h)Ylm{^ , (A.13) 

According to the first of \K.1\ . tfiis is equivalent to cos 02 being positive (negative) respectively. 
^"^^imiS) is the usual spherical harmonic function depending on the momentum angles. 
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is related to the radial wave function in coordinate space through 

Rniik)=ATT drr^Rni{r)ji{kr) , (A.14) 



where ji{kr) is a spherical Bessel function. The normalization is such that 

^3u _ roo 1.2, 



/^i5„,4*)r = 4.^°°i^ift.,wi^ = i . (A.15) 



We next turn to the explicit form of momentum wave functions for the ground state 
of the various atoms. For Hydrogen, the standard ground state wave function is 



with Zh = rriea being the inverse Bohr radius of the H-atom. In all cases involving 
integrals over the Hydrogen momentum wave function, we integrated over the range 
< ^ 3Zh, which describes very accurately the relevant momentum range of the 
bound electron. 

For the heavier atoms we follow ^T], based on the Roothann-Hartee-Fock (RHF) 
approach, also explained in PS.:- Accordingly, the radial momentum wave functions are 
written as 

Rnlik) = J2CjlnSjlik) , (A.17) 

j 

in terms of the RHF functions Sji{k) in momentum space, related to Sji{r) given in ref. 
[TT| through 

^ POO 

S,i{k) = A-K / drr^S,i{r)3i{kr) , (A. 18) 



in analogy to ()A.14|) . Since Sji{r) also depend on a parameter called riji, tabulated in 
the corresponding momentum RHF functions are: 



/ = 

n,o = l - 5,o(A:) = ^^^-^ , 

n,o = 2 ^ 5,0 A; = f ^ '° , 

~ , ^ 64Vl07iZ%\Z%-e] 
rijo = 3 ^ Sjoik) = - ' 
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~ 6AnZ%\5Zf,-lOZ],k^ + k^) 



Z = 1 

Uji = 2 ^ Sji{k) 



QA-nkZ]i^ 
64v^7rA;Zjf (5^2 - A;2) 

128x/T47rfcZ,^y^(35Zt, - 42^2^2 + 3^^) 



/ = 2 

= 3 ^ Sj2{k) 



128v^7rA;2zJf 
5(^^ + F)4 



1287rA;2#(7Z2 -fc2) 



The parameters {Cjin,nji, Zji) appearing in ()A.17l IA.19IIA.2T|) are given in Tables 
in [11]. Parameters {Cjin, riji) are dimensionless, while Zji are expressed in units of 
rrieOt = 3.73 keV in [TT] . 

We also note that, according to ()A.19IIA.2T| IA.17|1 and [TT], the highest Zjq for each 
atomic state determines the relevant range where the bound electron's momentum mostly 
lies^^. In most cases, this range is found to be < /c < 3 ma.x{ZjQ)mea. 
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Compare 



and the related variable Zh- 
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e-(P2) e-(P4) e-(P2) v{P-i) 

Figure 1: The Feynamn diagrams for z/ee~-scattering. 




Figure 2: The energy distribution of the final electron in z/gC" (a) or z/eC^ (b) scattering, 
at the rest frame of a free initial e~. Here Ee = is the final electron energy, and E^, 
[Ep) denote the incoming neutrino (antineutrino) energy Ei; see text. 
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Ee-me (eV) Ee-me (eV) 



Figure 3: The electron energy distribution in Vf, (a) or Vf, (b) ionization of if, and their 
respective ratio (c) and (d) to the corresponding distributions when the initial electron is 
assumed as free. Variables defined as in caption of Fig|21 
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Figure 4: The electron energy distribution in z/g (a) or z/g (b) ionization of He normalized 
to one e~ per unit volume, and their respective ratio (c) and (d) to the corresponding 
distributions when the initial electron is assumed as free; see Fig|21 Variables defined as 
in caption of FigO 
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Figure 5: The electron energy distribution in z/g (a) or z/g (b) ionization of A^e normalized 
to one e~ per unit volume, and their respective ratio (c) and (d) to the corresponding 
distributions when the initial electron is assumed as free; see Fig|21 Variables defined as 
in caption of FigO 
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Figure 6: The v^. (a) [z/g (b)] ionization integrated cross sections for H, He, Ne and Xe 
atoms divided by Z, and the cross section off free electrons at rest, as a function of the 
Vf. \pf. ] energy. The ratios of these ionization cross section, to the free electron ones are 
given in (c) and (d) respectively. 
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Figure 7: The energy distribution of the final electron in i/j^e' (a) or u^e' (b) scattering, 
at the rest frame of a free initial e~. Here Eg = E4 is the final electron energy, and E^, 
(Ep) denote the incoming neutrino (antineutrino) energy Ei. Identical results for 1/^. 
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Figure 8: The electron energy distribution in (a) or z/^ (b) ionization of A^e normalized 
to one e~ per unit volume, and their respective ratio (c) and (d) to the corresponding 
distributions when the initial electron is assumed as free; see Fig|71 Identical results for 
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